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Abstract 

Current and upcoming experiments should improve the sensitivity to Ea — >■ e^y decays by an order of magnitude. 
This paper assumes that one of the r e/37 decays is observed, and explores the structure and consequences of 
the required new flavoured couplings. In simple models (a low-scale seesaw, leptoquarks) it is shown that the 
^SJ . dipole vertex function is proportional to a product of flavoured matrices from the Lagrangian (a "Jarlskog-like" 

invariant), provided that the loop particles are weakly coupled to the Higgs. Secondly, if the dipole vertex function 
' has a hierarchical structure, this can imply that only some of the r — ^ e/37 modes can be observed, due to 

(n: the "approximate zero" implied by the bound on — !> e-y. The assumptions underlying this potential test of a 

O ■ hierarchical structure are discussed. 

^ j 1 Introduction 

Qh Measuring Cq — )■ epy would indicate New Physics near the TeV scale. Various studies of TeV New Physics models 

^ ' anticipate correlations in the Lepton Flavour Violating rates for Eq, — )■ e^7 and other processes. These predictions 
depend on the particle content and masses, the flavoured couplings of the model, and on the parameter space scan. 

■ This paper focuses on the flavour structure of New Physics: what can Ca — >■ e/37 tell us about patterns of lepton flavour 
, violation, with limited knowledge of the new particle content? 

■ After electroweak symmetry breaking, the dipole operator which induces ea — > 6^7, can be written, 
^" e / X 

a^■ -^{[XL]pc.ep^'"'PLeo.F^,, + [XR]0a.epa'"'PReo.F^,J . (1) 

' Lot: \ / 

where ct'"^ = 1(7'^, 7'^]/2. The dipole "vertex function"^ can be expressed as a dimensionful unflavoured constant, 

^ I ■ multiplied by a dimensionless "flavour tensor". A first question, is when and how does the ea — > 6^7 flavour tensor 
] reflect the flavoured couplings of the Lagrangian ("spurious")? We study, at one loop in two models — the type I 
• • ■ seesaw and leptoquarks — when the Ca — ^ 6^7 flavour tensor can be constructed, like a Jarlskog invariant [1], by 
contracting flavoured tensors present in the Lagrangian. 

Experimental searches for Lepton Flavour Violating (LFV) radiative decays impose significant bounds on the 
5— i ] fiavour changing dipole vertex function. The MEG [2] experiment constrains BR{fi 67) < 2.4 x 10~^^ and plans 
to reach a sensitivity of 10""'^'^ in the coming years. Babar and Belle have searched for r — >■ 6^7 decays, and impose 
BR{t ^7) < 4.4 X 10~^[3, 4] and BR{t ej) < 3.3 x 10^* [3] (see table 1). Super-B factories expect a sensitivity 
BR ~ 10~^ [5, 6, 7]. A second aim of this paper, is to explore whether Super-B factories, combined with MEG, 
can test the flavour structure of the dipole vertex function. So we make the non-trivial assumption that one of the 
T — e^7 rates is within the reach of the Super-B factories. This is a phenomenologically interesting scenario to 
envisage, because observing one r — > 6^7 rate gives two pieces of information, the second one being an "approximate 
zero" in the dipole flavour tensor, imposed by the /i — >■ 67 bound. We arge that, if the flavour tensor is hierarchical-'- 
and certain other assumptions are satisfled, then one of the remaining r — >■ 6^7 rates is suppressed below Super-B 
sensibilities. The assumptions required for this potential "test" of the hierarchical structure are discussed. 

Section 2 gives the Cq, — >■ 6^7 branching ratios, reviews various properties of the dipole vertex function X^^ij, 
and presents one-loop formulae for it. At the end of the section, is discussed the relation of X^^ij to the coefficient 
of the dimension six, electroweak gauge invariant, dipole operator. In the models considered in sections 3, which 
are a TeV-scale seesaw, and various scalar leptoquarks, the vertex function r turns out to be proportional to an 
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tThe ^L,R are arbitrary functions of coupling constants and masses from the Standard Model and beyond. They are refered to as 
vertex functions, to distinguish them from coefficients of operators in an expansion in 1/Apfp, where Apfp is a New Physics scale, 
-^meaning that it is dominated by its largest eigenvalue, analogously to [yu^Jlij — Vt^ti^^^i 
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BR 


current bound 


future 


T — >■ /i7 

r — >■ 67 


2.4 X 10-12 [2] 

2.5 X 10"^[3, 4] 
1.9 X lO-'^p] 


~ 10-13^ (MEG [2]) 

^ 10"^, (supcr-B factories[5]) 
~ 10"^, (super-B factories[5]) 



Table 1: Current bounds and hoped-for sensitivities to lepton flavour violating branching ratios, normalised to leptonic 
weak decays, as in eqn (2). 



"invariant" , when it is the coefficient of this dimension six operator. The definition of what we mean by "invariant" , 
can be found at the beginning of section 3. Section 4 reviews the argument that a hierarchical structure can be tested. 
Finally, the prospects of obtaining invariants and testing a hierarchy are discussed in section 5. A summary of useful 
loop results from a paper by Lavoura[8] appears in the Appendix. 

2 Notation and review 

For reviews of the expectations and prospects for lepton flavour violating processes, see e.g. [6, 7, 9, 10]. The 
theoretical predictions have been widely studied in popular models such as Left-Right models [11], multiple Higgs[12], 
supersymmetry [13, 14, 15, 16], Little Higgs models [17, 18] and extra dimensions [19, 20, 21]. The expectations of 
an A4 flavour symmetry, implemented in dimension six lepton flavour violating operators (so a New Physics model is 
not required), were explored in [15]. A more model-independent approach [13, 22], uses a "penguin-box" expansion 
[23] , to classify New Physics models according to the dominant effective interactions which they induce among ( three 
or four) SM particles. The various possibilities predict correlations [13, 16, 22, 24, 25] among observables (such as 
{g — 2)^ and n — )■ e7[26], eee and ^ e7[27], or /i — )■ 67 and /X6 conversion[28]). 

This paper studies flavour structure in the three — )• 6/37 decays, which are induced by the three-particle dipole 
vertex given in eqn (1). Electric and magnetic dipole; moments are neglected because they are flavour diagonal (see 
the expansion (3)), and other observables are usually neglected because they are induced by other operators. 

Lepton flavour violating radiative decays, Ca 6^7, proceed via the operator (1) at a rate given by 

r(ea e^Vcti'^) 25677* Gpm% 

where a, ^ e {6,/i,r}, and the whole paper uses the charged lepton mass eigenstate basis. Notice that in table 1, as 
in the remainder of this paper, the r — )■ 6^7 "branching ratios with tilde" are dcflned with respect to the leptonic 
r decay rate to facilitate comparaison with /U — > 67. Table 1 lists the current bounds, and hoped for sensitivities of 
running or planned experiments. 

The dipole vertex functions of eqn (1) satisfy Xjj = X/^: [X/jje^ induces jiji — > 6^7, [X^je^ induces /i/^ — > e/j7, 
[X/{]^e induces (/xl)"*" — >■ {eR)~^^, and [Xi]^e induces (hr)^ — (eL)+7. Since we are interested in — >■ 6^7 rates 
(and not CP violation), we only need to consider the a> j3 components (lower triangle) of [XijJ^a and [XlJ^q,, which 
are independent. The dipole vertex functions can be written as : 

Xfl = t; (ci?Y^ + SlY^ + Ti^ ) ^ (3) 

V 2 
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where v = VIA GeV is the Higgs vev. The normalisation by 2-\/2Gf is concrete, however, the masses of the new 
particles which mediate Ba 6^7 will be assumed to exceed the weak scale v. Ye is the (diagonal) charged lepton 
Yukawa matrix, cl and cr are unknown, unflavoured constants, and S^, S/j and T^^r are new dimensionless matrices 
in flavour space, which we will refer to as "flavour tensors". Section 3 studies when they can be constructed by 
multiplying together flavoured matrices that one finds in the Lagrangian. Si and are hermitian, respectively 
induced by New Physics which couples to doublet leptons, or to singlet leptons. Tr = is an arbitrary matrix (as 
was X), Tr is defined like the Yukawas to have doublet-singlet index order, induced by New Physics interacting with 
doublet and singlet leptons. The 'S.l^^r are labelled by the chirality of the incoming charged lepton; the Sl.r have 
the opposite chiral label, because the incoming chirality is fiipped by a charged lepton Yukawa coupling. 
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Figure 1: The first blob corresponds to [Si]e^, the second to [Sijje;:^. On the second row, the blobs correspond to 
[Tfl]e^ and [TlI^. 

The decomposition of eqn (3) corresponds to flavour diagonal + flavour changing new physics that is not chirality- 
flip (the Sl.r terms therefore contains Vg) + flavour changing chirality flip new physics represented by Tl,r- The 
first term oc Ye, whose coefficient can describe the magnetic and electric dipolc moments, can not induce the flavour 
off-diagonal — >■ e/37. This term is nonetheless included, so that the diagonal elements of our new flavour tensors 
are not required to parametrise flavour diagonal new physics. The lower triangles of Sl.r, T^l and correspond 
respectively to the four diagrams of figure 1. We have neglected the possible Y^S^ contribution to Xj;,, and the YgSjj 
contribution to X/j, because they are proportional to the outgoing charged lepton mass m^. 

In the presence of either the S^.ij or T flavour tensors the ratios of branching ratios of eqn (2) become (recall 
that 4Gf/V2 = 1/^2, for v = 174 GeV): 

BR{l^^lp^) = —{\Sl0c.? + \Srp^\^) 

~ 1.4 X 10-2 + |Sk;3„|2) (4) 

BRQ^^l^j) ~ l^{\TL^^f + \TRp^f) . (5) 

If the the dimensionless couplings in Txfia and Sxfia are taken ^ 1, such that Sxfia ~ ^"^ I^npfg ^"^^ Txpa ~ 
/ A?j^pp(j, then Super-B factories could probe to a new flavour changing scale A^pfc ^ 10 TeV for Sl,^ and 
AjvpFC ^ 60 TeV for T. 

Convenient formulae for the [Xx]/3a vertex function, generated in one-loop diagrams involving a boson and a 
fcrmion / (see e.g. flgures 2 and 3), are presented in the paper of Lavoura [8]. For a scalar S with Yukawa-type 
interactions 

^(Lavoura) 

= SfiXLaPL + \RaPR)ea + h.C. (6) 

where ea- is a charged lepton of flavour a, in four component notation, Lavoura obtains 

^ = -i^^[Qsm+Qsm] (7) 
^ = -^^[Qfm+Qsm] (8) 

^ = + (9) 



(10) 



where the electric charges satisfy Qg — Qf = 1, the dimensionless functions k,k,kf and kf of the mass ratio t = mj/mg 
are given in the Appendix, and 

W/3a = A^^Aia , P/3a = A^^A/Jq. , ^^q. = X*LpXRa , C/3a — X'^i^Xlu ■ (H) 



§if both are present, there are ST interference terms 



3 



Lavoura also gives formulae for a vector and fermion loop. Suppose that W bosons interact with clo- and a neutral 
fermion / via the Lagrangian 

£,ec = A'^^WjYeLa + h.c. + SCf'h (12) 

where Aj^^ = 5/-\/2xmixing angles, and 6Cf'h is the additional interactions, in Feynman 'tHooft gauge, of the Higgs 
goldstone components. Then Lavoura gives 

= -^Qwm , where p/sa = A^^A^^ (13) 

y is given in the Appendix, and [SL]i3a = [Tnl/Sa = [TL]i3a = 0. 

In the SU{2) X U{1) invariant Lagrangian, there is a dimension six operator O^^ , which contributes to the photon 
dipole operator, and which is a linear combination of hypercliarge and weak SU(2) operators: 

OZ = {l0H)a''-'Pne^F^, = cos dw CWu^" PrBo^B^, - sm0w{i()T^H)a'''' Pne^W^, (14) 

where £ is an SU(2) doublet, cr'"' = i[7'*,7'^]/2, and {r*} are the Pauli matrices. If the operator of eqn (14) appears 

in the Lagrangian normalised as £ D e2\/2GFC^^2o^a + h.c, then C^^^ contributes to [X^J^a (see eqn (3)). Notice 
that there could also be higher dimensional dipole operators, such as 

O'^f' = {H^H)(i0H)a<'-PRe^F^, (15) 

with additional Higgs legs and suppressed by more powers of the flavoiu- changing New Physics scale Ajvp_fc- There- 
fore, C^^) corresponds to the 0{1/ h^f^ppo) terms in X/j, and the coefficients of O'^'^^ include C'(to|m/-^jvpfc) 
contributions^ (where msM is a SM mass). In this paper, we assume that there is a New Physics contribution to C^^^ 
and we will be interested in identifying cases where C'^' ~ Xjj. 



3 Invariants 

For the purposes of this paper, invariants are tensors in flavour space, often a scalar or a matrix, obtained by multiplying 
flavoured matrices from the Lagrangian ("spurious", in Minimal Flavour Violation [29] langauge). They are interesting 
for Beyond the Standard Model physics, because, if a sufficient number of invariants are measured, the flavoiircd 
matrices of the Lagrangian can be "reconstructed" by simple matrix multiplication. Like the original invariant of 
Jarlskog [1], these invariants are also an elegant way to avoid dependance on unphysical or unknown basis choices in 
the Lagrangian. So the invariants in this paper are allowed to have flavour indices, but these must be in the mass 
eigenstate basis of known particles. 

It is clear that that the formulae (7) - (10) and (13) have the makings of "Jarlskog-like invariants", at least in the 
case where the /c-functions can be approximated by a single term. For instance, if rn/ <C mg, then Sl.r ~ Am^^A'^. 
This would be a coefficient of the dimension six operator given in eqn (14). The remaining terms from the k functions, 
DC {mf/ms)^ for n > 1, would give coefficients for dimension > 6 operators which give more suppressed contributions 
to the dipole vertex functions. The first aim of this paper, is to identify when the flavour tensors Sl.p, and Tr 
are invariants. Various questions come to mind about this connection: what happens to logarithms? What if both 
particles in the loop are flavoured? What happens when there are several new particles with slightly difi'erent masses? 
Whether the vertex functions are invariants is essentially a "top-down" question. However, it is phenomenologically 
interesting, because we would like to reconstruct the fundamental Lagrangian from observables, so it is useful to know 
when the connection is simple. 

3.1 The non-supersymmetric Seesaw 

The type I seesaw [30], (without supersymmetry), is simple model in which to study the invariant question, because 
it only contributes to S^, and because only one of the particles exchanged in the loop (the neutral fermion) carries a 
generation index. The singlet neutrino masses are assumed to be at the TeV-scale[31], as can for instance arise in the 
"inverse seesaw" [32]. Some recent constraints on such models can be found in [33, 34]. This mass scale is chosen so 
as to obtain detectable Bq, — > e^7 rates [35, 36], while maintaining the new mass scale above the Higgs vev v. (The 
seesaw behaviour works for singlet masses ^ TeV [37], but our discussion will not apply.) 

^The coefficient of O"''^^ would also include the contributions oc to the mass^ of the new particles. However, if we propagate mass 
eigenstate new particles, such contributions would remain resummed on the denominator. 
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Section 3.1.1 approximately diagonalises the neutral mass matrix, and checks that the known epj amphtudes 

are reproduced " using the formulae of [8] for the neutral fermion and W loop diagram. Section 3.1.2 studies the 
conditions under which the amplitudes are proportional to invariants. 



3.1.1 The model 

Consider adding r heavy {Mj ^ TeV) singlet fermions N , with a majorana mass matrix M, to the three generation 
SM. In a simple inverse seesaw model, there would be r = 6 singlets, arranged in three pairs of opposite-sign, but 
almost-equal magnitude, masses. The leptonic Lagrangian, in the mass basis of charged leptons and heavy singlets, 
will be: 



£ = KdocKNi-^NiNi-yJ„H*^e„ + h.c. 



(16) 



where = e(i?„ )* — u{H^)* , and two Higgs doublets are allowed for. The neutral mass terms will be 

Cmass ^ -Xc,iV^{H^*)Ni-^NiNi + h.C. 

and the resulting (3 + r) x (3 + r) majorana mass matrix can be diagonalised with a unitary matrix X: 



X 







XVu 

M 



X' 



Dm 

Dm 



where £)„ 



-vJ.U^XM-^X^U* 



(17) 



where Dm = diag{Mi, ..Mr}, and we suppose that it is sufficient to obtain X to order because O^'' of eqn (14) 

is a dimension six operator. Assuming that det [M] ^ 0, because all the singlets are heavy, and taking U to be the 
usual leptonic mixing matrix, gives 



X = 



Jjt 
1 



1-dl 



-AM-i 
I -62 



where the 6i ensure that X is unitary at 0(1/M^): 

Si = ^A[MtM]-iAt 



62 = ^[Mt]-iAUM-i 



(18) 



(19) 



A flavour eigenstate neutrino, participating in a vertex with the W and a charged lepton, can therefore be expanded 
on light (vi) and heavy (nj) four-component majorana mass eigenstates as 



(20) 



3.1.2 To get invariants? 

The contribution to the dipolc flavour tensors is obtained by summing the diagrams where a W boson and a light or 
heavy neutral fermion is exchanged. The relevant diagrams are given in figure 2. We neglect the light neutrino mass'^ 
contribution because it is parametrically of 0(A**), to be compared to the O(A^) contributions from heavy neutral 
exchange and from the "non-unitarity" of the light neutrals. 
Combining eqn (13) with eqn (20) gives 

l^[X[M^M]-'xnp^{-y,{0)+yi{Mj/m'w)) 

-^[X[M''M]-^X'']fta + - heavy + light in loop (21) 

where Vu/v = sp and y is given in eqn (49), and only the 0(1/M^) terms were retained. This is consistent with 
calculating the diagonalisation matrix X to 0{1/M'^). 

" Despite comments in [34] , the correct amplitudes can also be obtained by first matching the unbroken gauge theory at the heavy scale 
M onto the dimension six "kinetic" operator iH pH*i of [20] and the dipole operator of eqn (14), and then matching at mw, in the broken 
theory, onto the dipole, including the kinetic operator in the loop. At the high-scale matching, it seems that the not-lPI diagrams with a 
Higgs insertion on the external leg should be included. 



y2 
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Figure 2: One- loop diagrams that contribute to the dipole vertex function in the seesaw model with broken electroweak 
symmetry. H is the goldstone. A Higgs leg, which is not drawn, can attach to the /i or m line. 



As hoped for, equation (21) is an invariant constructed from flavoured matrices from the Lagrangian. This is the 
case despite electroweak symmetry breaking and the separate contributions from the scales M and m^y. However, 
there are 0{v'^ /M^) contributions from higher order terms in the expansion of y. From eqn (49), one sees that heavy 
and light neutral exchange will give contributions with coefficients proportional to 

t;2[A[MtMMtM]-iAt] ... , [XM'^ \^ \* M'^* X'^] , v^[X[M^ M]-^X'^XM-^X^X*M-^*X'^] ... (22) 

Since the eigenvalues and eigenvectors of the neutral mass matrix were only obtained to 0(1/M^), the numerical 
factors multiplying the 0{v'^ /M^) terms arc unknown. The coefficients are still invariants, which seems reasonable 
when the flavoured particles in the loop all have masses beyond current experimental reach. Fortunately, there are no 
logarithms at dimension six. 



3.2 Scaleir leptoquetrks 

The second toy model is various scalar leptoquarks. For updated low energy bounds, and references to earlier works, 
see e.g. [38]. The current collider lower bound on the mass of a leptoquark decaying to first generation leptons and 
quarks is 660 GeV [45], based on 1 fh~^ of data. See [44] for earlier collider mass bounds. 

The scalar leptoquark contribution to the dipole vertex function is more complicated than the seesaw model, because 
they can contribute to Sl,^r, or T (because they can interact with both singlet and doublet charged leptons), and 
in that both they and the quark in the loop can be flavoured. 

Subsection 3.2.1 gives the Lagrangian for baryon and lepton number conserving scalar leptoquarks. and the contri- 
butions to Cq, — >■ e^7 of two SU(2) singlet leptoquarks upon which we focus. These are representative of the possibilities 
(see table 2), because the Sq leptoquark can interact with singlet and doublet fermions. Then subsection 3.2.2 discusses 
the prospects for obtaining invariants. 



3.2.1 The models 



Consider scalar leptoquarks, with renormalisable B and L conserving interactions, which can be SU(2) singlets, 
doublets or triplets. In the notation of Buchmuller,Riickl and Wyler[39], these can be added to the SM Lagrangian 
as: 



SoiX-LSoiinq" + XrsoPu") + 

+ (Als2^" + ARS2eg[iT2])52 

+ (Xi^sJlT2Tq'')S3 + h.c. (23) 

where the As are 3x3 matrices with index order lepton-quark, and {tj} are Pauli matrices, so iT2 provides the anti- 
symmetric SU(2) contraction. In this Lagrangian, the leptoquark leaves the vertex into which enter the leptons. This 
is converse to Lavoura conventions, where scalar and lepton both enter, and the internal fermion leaves. Comparing 
eqn (23) with eqn (6) gives the parameters listed in table 2. We will also allow the possibility of three generations [40] 
of leptoquarks, in whiich case / : 1..3, and Xj^ is a three-index tensor. 

The SU(2) singlet leptoquark Sq, which couples to d'-'e, contributes to via the two diagrams summarised by 
the left diagram in figure 3 (the photon can attach to S or f), where the internal fermion / can be a singlet or 
6"^ quark. For a 5*0 leptoquark, eqn (7) implies that = T^ = Tr = 0, and 



[S 



R\l3a 



E 



^(M-)+4Mx)) 



1 



d^RSo^R^M 



leptoquark without generations (24) 
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Si 


1/3 


-2/3 


[ALSo]/9g[ALSo 


* 
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[AflSo 


PI [AkSo 


*af 


[ALSo]/3/[A_R5o]a/ 


[A_RSo]/3/[ALSo]a/ 


Sl {lower) 


5/3 


2/3(u) 


[AlS2]/3/[Ai,S2 


«/ 


[AflSa 


/3/[AiiS2 




-[ALS2]/3/[AilS2]*/ 




Sliupper) 


2/3 


-l/3(d) 







[AijS2 


/3/[AiJS2 










sl {lower) 
{T^comp) 


2/3 


-l/3(dij) 


[AlSsI/S/I^LSs 


a/ 















1/3 


-2/3 {u") 


[Als3]/3/[Az,s;3 


a/ 
















Table 2: Parameters for obtaining the Ca — >■ e^7 amplitude, induced by the scalar leptoquarks of the left colomn, 
using eqns (7) - (10)). a and /3 are lepton flavour indices, and / is the internal quark index. See flgure 3. 




Figure 3: One loop diagrams contributing to — )■ e/37, mediated by a scalar S, and a fermion /. In the diagrams on 
the left, the scalar S interacts either with doublet or singlet charged leptons, so the Higgs insertion is on an external 
leg. In the diagrams on the right, S interacts with singlet and doublet charged leptons, so the Higgs insertion can be 
on the internal lines. 



where x = m^/M"^ <^ 1. Inside the square brackets after the first equality, psa is from table 2 but not summed over 
quark flavour q. Assuming x is negligeable, k ^ 1/6 and k 1/12, which gives the second equivalence; this means 
that the coeflacients oyimension eight operators, in this model, axe irrelevant. 
In the case where So has generation indices /, one obtains 

- -j5I[^ko7v^^flkl''" = "i[^fl^o-OM'A^]/3o (25) 
/ s,i 

where after the first equality, there is a [A^] matrix with indices in lepton and quark fiavour spaces for each fiavour / 
of leptoquark, and after the second equivalence, the leptoquark index is also implicit {Dj^f = diag {M| ^--M^ ^} is the 
diagonal leptoquark mass matrix). A^g is in the charged lepton, down quark and leptoquark mass eigenstate bases, 
and the quark flavour sum has again been done neglecting quark mass effects oc mg/Mj ^. 

Consider now the singlet Sq coupling to r^l.r, without any flavour index of its own, because it is straight- 
forward to add a flavour sum on the leptoquarks to the formulae below, couples to up-type quarks, so^ = u, c, t, 
and the 0{m\, m\)/M'g terms in k and k can be neglected, as were the quark mass terms in the case of Sq. We keep 
the top quark mass dependance. For x = mj/Mg, t = m^/Mg, one obtains [8] 

^ = -E^HM.) + fc(.)) 



[AijSolfAijSo]^ 



/3a [^RSoht[^RSo]lt ft^-M^ + 9t t{t + 2)\nt 



8M^ 2M| V 4(t-l)3 2{t-l)^ 



(26) 



s 



f 

/3a [XLSohtiXLSoVat ft^ -St^ + t{t + 2)lnt 



[[Xls„][XlSoV 



8M§ 2M| V 4(t-l)3 2(t-l) 



(27) 
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2M|w V2(i-1)2 (^-1)^ 



[AflSo]^/m/[Az,sJ*^. (7 ii^f 
~ 2Mlv M| 



(28) 
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where p, co, ^ and C are from table 2 with no sum over q. 



(29) 



3.2.2 Getting invariants 

For the So, which interacts with charged leptons and down- type quarks, the dipole vertex function is given by eqn (24), 
or (25), respectively for the cases that So does not, or does, have generation indices. It is an invariant, constructed 
from spurions from the Lagrangian: 

y2 = ~'^l^RSo^Ai^RSo^^°' ' leptoquark ^0 with generations (30) 

where the quark flavour sum is implicit. Since the down-type quark masses are small compared to M~, dimension 
eight and higher operators are irrelevant, the dipole vertex function in this case is unambiguously an invariant. 

The dipole vertex functions generated by the Icptoquark 5o, which couples to up- type quarks, are not so obviously 
invariants due to the large top mass. For Sl.r, if ^■^i3tXat''^t /^s ^ [''^^^]/3a, the second term in eqns (26) and (27) 
can be ignored, and the Sx are approximately the first term, which is an invariant: 



[AxSo][^XSo]^ 



8M| 



+ ... (31) 



However, if the terms of order rn^/Mg should be included, they correspond to dimension eight operators with invariant 
coefiicient a Ay„yjA'^ (in the one-leptoquark-generation case). This ressembles the seesaw case, where the coefficients 
of dimension eight operators were potentially dangerous, because not nccccssarily suppressed by small couplings. 

In the presence of Al and Xr, a contribution to T can arise, propertional to the internal fermion mass, giving the 
lowest order coeflacient proportional to invariants of the form XrY^X''^ or AlF„aJj. So one obtains 



-L\0a 



= -(^[Afl5o(cii;^ + Fjln(rjK))A^w) +... , 
^ = -(^J>^LSo{ciYu + Yuln{YuY;i))X^j^^Jpo) + .■■ , So with Xl,Xr 



(32) 



where Ci = 7/4 + ln(v2/M|), and the quark flavour sum is implicit. The logarithm is unfortunate, from the perspective 
of reconstructing the new flavour structures, because one needs its coefficient to reconstruct the new fiavoured matrices. 
For instance, imagine that [Ai,]/3q was determined in some other process, and that the pattern of Cq, — ?• e^7 decays 
indicates that they are mediated by a hierachical flavour tensor with chirality flip on the internal line. Then to obtain 
Xr from Ca — >■ 6/37 data requires knowing the coefficient of the logarithm. Two positive feature of logarithms are 
that one can often guess their presence, and, as discussed in section 3.2.3.1 of [10], if one knows the coefficient, the 
logarithm does not obstrct the reconstruction of new flavour structures. 

Eqn (32), like eqn (31), neglects the dimension eight terms of order XrYuY^YuXlv^ /Mg and XlY^Y^YuXrv'^ /Mg, 
which are not small for the top Yukawa. 

Finally, it is clear that the dipole coefficient induced by will only be proportional to a single invariant, if one of 
^L.R, Tl or Tr dominates over the others. See the comparaison of their relative sizes in section 4.3.2. 
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4 Testing hierarchies? 



It is well-known, that if the flavour tensors Sl.r or T are dominated by their largest eigenvalue, then they are 

compartively predictive. The idea is that, if one of the r — > e^7 decays is observed, then the current bound on 
/Lt — >■ 67 implies a small mixing angle in the diagonalisation of the flavour tensor, and this small angle suppresses 
the other r — >■ £'j decay. This section is somewhat independent of the previous section, because the flavour tensors 
contributing to the vertex function can be hierarchical without being invariants. Nonetheless, subsection 4.3 discusses 
when hierarchical invariants are obtained, because ideally, one wants to reconstruct the spurions in the Lagrangian 
from measurements of the vertex function. 



4.1 Defining iiiercirciiies 

We define a matrix, in particular the flavour tensors S^jj, T^, and T^^, to be "hierarchical", when the off-diagonal 
elements are dominated by the largest eigenvalue. A hermitian matrix S (for instance Sr or SL),can be written 

S = V^DsV , = diag{si,S2,S3} , Si < < S3 (33) 

and S3 will give the largest contribution to Sj^a, when 

n>l^la , 1^2>S2^^2a < ^^3>S3^3a (34) 

Eqn (34) is satisfied for eigenvalues Si,S2 and S3, if the mixing angles are large enough **: 




This relates elements of a colomn in V: the bottom row [Vj^a has to be bigger than a fraction y/sl of the elements in 
the colomn above it. 

The flavour tensor = T^, like a Yukawa matrix, can be diagonalised with independent unitary matrices on the 
left and right: 

Tn = vlDTVR , Dt = diag{ti,t2,t3} , h<t2<ts (36) 

If the matrix elements of Vl and Vr satisfy eqn (35) (replacing Sj U), then, just as in the case of Sl.r, the 
off-diagonal elements [T^J^q,, and [Tfl]^c will be dominated by t^. 



4.2 Observing a hiercirchy 

Suppose initially that New Physics induces only one of Sl or Sr (the possible presence of Si, and Sr, or Tl,r will 
be discussed later). Then we start with two assumptions 

1) the fiavour tensor Sx is hierarchical 

2) one of the t ^ ep-f decays is observed at Super-B factories. This is neccessary for them to be able to test anything. 
To be concrete, suppose that r — >■ /X7 is observed. 

Combined with table 1, these assumptions fix 

3 X 10-3 ^ \V3*^ssV3^,\ ^ 10-3 (37) 

where S3 ^ A^w^/A'/^^^^ < w^/A^^SAf (assuming that Sl.r is induced at dimension six by New Particles with couplings 
A < 1). In addition, the non-observation of /x — t- 67 puts an upper bound 



BRifi^e^) _ l^e^P _ 1^3 



3e| 



|2 



^10-^ . (38) 



With this bound on Vse, the only way that r — )■ 67 could detectable rate at Super-B factories, is if — V^e are both 
small. However, from eqn (37), this would require S3 > .1, or Mrsm ^ 500 GeV, which is a low scale for new particles 
with 0{1) couplings to charged leptons^^'. If we make the further assumption 



M 



3) that S3 ~ —2 <C .1 , that is, Mbsm ^ 500 GeV 



BSM 



** Equation (35) implies (34), but since eqn (34) is a lower bound on a product of elements of V, it can be satisfied when one of the Vij 
does not satisfy eqn (35). 

ttThe 7 TeV LHC, with the current 5 fh~^ of data, is surely sensitive to leptoquarks of mass ^ 600 GeV — even if they decay to rt [41] 
— and possibly also to heavy singlet fermions 
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then we can conclude that Super-B factories can falsify the assumption of hierarchical couplings in or S^, by 
detecting both r — > fi^ and r — ej. 

An alternative way to see the fj, ^ bound, for a hierarchical flavour tensor, is illustrated in figure 4. The 
paremeter space of the three elements iVsap, who satisfy l^^l^ ~ 1' ^^"^ equilateral triangle. The blue 

hyperbola corresponds to T{t — ^7) cx |V3^V3t-P visible at a Super-B Factory. Then the green plane corresponds to 
the bound of eqn (38); allowed points live on the intersection of the red and blue surfaces, and on the far side of the 
green plane. The hyperbola corresponding to a fixed rate for r — >■ 67 is not drawn (it would be in the vertical and 
I Vsep plane, intersecting the blue hyperbola). The green and red planes are independent of the (unknown) eigenvalue 
S3, which is oc /M'^gj^^. The hyperbolae do depend on S3, and move away from the axes as S3 decreases (or as the 
branching ratio increases for fixed S3). For the value of S3 chosen in this plot, it is clear that the t — )• 67 hyperbola 
corresponding to a rate similar to t /X7, would live on the near side of the green plane, so a Super-B Factory cannot 
see r — >■ 67 and t — >■ /X7. 




Figure 4: A graphical representation of the claim that a hierarchical flavour tensor (see section 4.1) predicts that 
a Super-B factory should not see r — 67 and r — /17. The axes are {IVsep, jVs^p, IVstP}, which satisfy jVsep -|- 
|^3/iP + |^3tP = 1 represented by the red triangle. The vertical axis is iVsT-p. The blue hyperbola corresponds to 
r(T — )• /i7) cx iVg^VsrP visible at a Super-B Factory. The green plane corresponds to the /x — )■ 67 bound of eqn (38). 
Allowed points live on the intersection of the red and blue surfaces, and on the far side of the green plane. A hyperbola 
(not drawn) corresponding to r — >■ 67 at a rate similar to r — >■ /i7 would live on the near side of the green plane, see 
discussion after assumption 3), so a Super-B Factory should not see r — >■ 67 and r — ?> fi'y. 

The assumption that only one of Sl or Sr is present is minimal, and perhaps justified by the absence to date of 
new TeV-scale particles. However, it can in principle be tested: if the decaying lepton is polarised, then by angular 
momentum conservation, the angular distribution of the final lepton will depend on its chirality [9]. Therefore, if 
r — >■ /X7 and r — > 67 are both observed at Super-B factories (where polarised ts may be possible[6, 7]), one can hope 
to determine, from the final state angular distributions, whether they are induced by Sl and/or S^j. If the final state 
angular distributions are different, then this is compatible with hierarchical Sl and S^j; if both r — > /i7 and r 67 
have the same final state angular distributions, then this is incompatible with hierarchical couplings in S^.t?. So one 
more condition is required to test a hierarchy in the 8^ ^: 

4) in decays of polarised rs, the angular distribution of final state leptons can be measured 

However, New Physics which induces and S^; may also induce Tx, as was the case for the leptoquark Sq. 
Consider a hierarchical T, neglecting at first the S^.^. The Cq, — ^ e/37 branching ratios will therefore be parametrised 
by t^, and the third rows of Vl and Vr. This makes five parameters (as opposed to three for Sl or S;?), so the polarised 
branching ratios are required to hope to test the hierarchical pattern in T. Suppose, for instance, that tl — > /ii?7 is 
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X 
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TL 6^7 
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X 
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Table 3: If the polarised decay indicated by ^/ is observed, then a hierarchical pattern in the operators listed in the 
top row predicts that the decays indicated by x should not be seen at Super-B factories. The blank spaces mean that 
no prediction is made. This table depends on the assumptions 1-4 listed in the text, see in particular assumption 3). 
The operators Sj, and can induce either of their two colomns, T can induce any two non-conflicting colomns. See 
the discussion at the end of section 4.2. 



measured. Both BR{iil e-Rl) and BR{fj,R — >■ 6^7) are small, so, in particular 

BR{rL^mi) M I^lV^hsmP ^ ' 

so with assumption 3) above, BR{tr c-lj) \Vl^^VR3r\'^ should be beyond the reach of Super-B factories. Similarly, 
if tr^ Cl^ was measured, a hierarchical structure in T predicts that tl — ^ fiR'j is inaccessible. 

Table 3 summarises which processes should not be seen, given a polarised decay induced by a particular flavour 
tensor. For instance, if one r — > ex/sj decay is observed, this is consistent with a hierarchical Sx or T. If both one 
of the Tr — >■ 6^/37, and one of the tl — >■ cr^j are observed, this is consistent with simultaneous hierarchical Sr and 
Sl- If the particular combinations, tr ^lJ and tr clJ, or tl — >■ ^lRJ and tl — >■ gr'j, are observed, this is 
consistent with a hierarchy in T. So hierarchies in the various flavour tensors predict that various rates are suppressed. 
However, if all the r — >■ e/37 decays are observed in all polarisations, this is consistent with the simultaneous presence 
of hierarchical Sl.r and T. That is, a hierarchical pattern in T or in the Sl.r make opposite predictions: if, for 
instance, tr — > jiLj is observed, this is compatible with S^, in which a hierarchy imposes that tr — > 6^7 must not 
be seen. However, tr ^ ^l1 could also be induced by T, where a hierarchy is consistent with tr — >■ cl^, but not 
Tl eRj- Therefore, with the assumptions made so far, the hierarchical structure cannot be falsified. In addition, 
one must assume that 

5) either the Sl.r, or T is dominant. 

Then it is clear from the table that a hierarchical structure can be falsified: for instance, if tl — >■ /Ul7 and tl — cl^ 
are observed, this is inconsistent with a hierarchy in Si, (or Sr). 

4.3 Hierctrchies and invciriants 

This subsection makes contact between flavour tensors as invariants, and the prospects of testing their hierarchical 
structure. Section 4.3.1 discusses the conditions which the spurious should satisfy, such as to obtain a hierarchical 
invariant. Section 4.3.2 uses invariants to study whether the leptoquark So, which can induce Sl.r and T, satisfles 
assumption flve for testing a hierarchy. 

4.3.1 Getting a hierarchy in Sl.r 

Ideally, the hierarchy condition eqn (35) could be translated to conditions which the spurions should satisfy to ensure 
a hierarchical invariant. However, this phenomcnological approach amounts to expressing elegant invariants in explicit 
basis-dependent form, which defeats the purpose of invariants, and risks to produce opaque interlinked bounds on 
unknown eigenvalues and mixing angles. Instead, it is simple to study the issue in the wrong direction, that is, start 
from hierarchies in the spurions, and enquire when then are transmitted to the flavour tensors. 
Two simple cases are: 

1. degenerate eigenvalues Mj of the matrix M of the New Particles 

Consider to be concrete the seesaw case. If A is written, in the charged lepton and heavy singlet mass eigenstate 
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bases, as V^DxVr, then for degenerate singlet mass -squared |M/p, the hierarchy condition on S eqn (35) is 
trivially applied to A''' A in the charged lepton mass eigenstate basis: 

[VL]j0>^[VLy i<j (40) 

To explore the area around the degenerate mass limit, one can write 

[XD-^^X'^W = ^ (^[AAt],. - [A]ff J ^ A:,j (41) 

and assume that the first term satisfies eqn (35). Then the second term is negligeable if 

Mj - Ml < (42) 



2. flavour- independent vertices 

If, in general, A = Vj^D\Vr, then for Dx proportional to the identity matrix, the hierarchy condition of eqn (35) 
applies to V = Vj^Vr and D^. 

Illuminating formulae for the case where there are hierarchies in A and M are not obvious to find. 

For 5*0, without generations indices, S/j will have hierarchical form if A|j.A-^ = V,,D'jy^ is hierarchical on its lepton 

indices (equivalently, if Vg and satisfy conditions (35)). If Sq has generation indices, then for three leptoquarks of 
mass squared Mf, the same conditions as discussed for the seesaw are sufficient to obtain a hierarchical flavour tensor. 
Similar comments would apply to the leptoquark So with couplings Xl or A^. 



4.3.2 T versus Sl.r for the 5*0 leptoquark 

The fifth condition for testing/confirming a hierarchy, is that only one of T or Sl.r makes a significant contribution 
to T — >■ epj. The leptoquark Sq, with both couplings A/, and Xr, induces both structures of matrix element, so it is 
interesting to compare the relative size of T and r. 

First, notice that constraints from two-quark- two- lepton contact interactions apply directly to the magnitude of 
terms contributing to the invariant. This is a useful feature of invariants. The bounds from r — > epn [38] impose 
that the u quark contribution to [Sl.r]i3t is less than 2 x 10"'', so insufficient to generate an observable r — >■ e^7 
rate (see eqn (37)). On the other hand, a diagram with c quarks in the loop can induce an observable r — >■ e/37 via 
Sl.r or T; there are 0(1) bounds on [Sij]/^,- fron Z — >■ e^r, contraints from measuring Vcs on [Sl]/3t which are an 
order of magnitude smaller, and finally bounds from Ds — >■ TUp on the c-loop contribution to T which (just) allow it 
to contribute an observable r — e/37 rate. Finally, it is clear that loops involving a top quark can contribute to an 
observable rate via at most two of Si r and T. This is because one cannot simultaneously obtain [Ai,]*([AL]^tm^ ^ 
[XR]Tt{^R]^tfnT ~ [XLjTti^Rlfit'nit ^ [XR]*f.[XL]i3t'mf So unfortunately, So leptquarks interacting principally with a c 
quark could have hierarchical couplings that generate similar rates for the four processes r e^^R-y and r — >■ Hl,r1- 
This is illustrated in figure 5, which shows flat asymmetries (solid) and Aclr (dotted), defined as 

A, - ^ ~ ^ Ac^^ - ^ - ^ gfi^) (43) 

BR{t M7) + BR{t 67) ' BR{t 67) 

An asymmetry for muons A^lr can also be defined; it has the same distribution as Aclr- The histograms are 
obtained by varying the six couplings Xna Xrcc (for £ € {r, /i, 6}) between .01 and 1 with a log prior, and computing 
the asymmetries for all the points which give a total rate for r 6^7 below the current bounds and within the reach 
of SuperB Factories. The normalisation of the vertical axis is arbitrary. 



5 Discussion 

The observed neutrino masses demonstrate the presence of New Physics in the lepton sector. However, we do not know 
what it is. One approach to this problem is to calculate in motivated models, thereby restricting their parameters and 
identifyng their predictions. Alternatively, one can try to "reconstruct" the New Physics from data. This paper aims 
to follow the second, more phenomenological philosophy. 
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Figure 5: This plot illustrates the importance of assumption 5) for testing hierarchies. An So leptoquark coupled 
principally to c quarks would induce hierarchical flavour tensors, but can give similar rates to the four final states 
eL7, efl7 /xl7, and i^r'J. The sohd and dotted lines are respectively the asymmetries and Ablr (see eqn (43)), 
obtained after scanning over allowed couplings (see text around eqn (43)). 

5.1 Invariants 

If the masses and couplings constants of the new particles were known, the dipole "vertex functions" ^l.r of eqn 
(1), which arc observable, could be calculated to arbitrary accuracy. However, since wc know neither the new particle 
identities, nor their mass scale, nor the structure of their flavoured couplings, it would be convenient to be able to 
separate these miknowns in ^l.r- So the first aim of this paper was to study, for a variety of new particle contents, 
whether it was possible to fac;torizc ^l.r as the product of a dimensionful, unflavoured constant, multiplying an 
"invariant" . We refer to flavoured coupling matrices from the Lagrangian as "spurious" , and define invariants to be 
products of spurious. Since X^^^ is a dimensionful matrix in flavour spaces, it can always be written as a dimensionful 
factor multiplying a dimensionless tensor in flavour spaces, which we call a "flavour tensor" (see eqn (3). The crucial 
question is whether this flavour tensor is an invariant. This is an interesting question for phenomenology, because it 
is comparatively simple to obtain information about spurious (the Lagrangian flavour structures) from invariants. 

In this paper, we only consider one loop coontributions in New Physics models without a conserved parity. This 
means that in the loop, there is always a Standard Model particle, and a new particle. In such models, it seems that the 
question about invariants reduces to whether the 'Kl^r are dominated by the coefficient of the dimension six operator 
of eqn (14). This claim is based on comparing the formulae for the X^.a, to the Effective Field Theory expansion in 
1/Anpfc, where Ajvpfc is the scale of new flavour changing physics in the lepton sector. In eqn (3), the ^l,r, of mass 
dimension -1, are written as a sum of dimensionful constants ~ AG pv / {IQn-^ y^) , multiplying dimensionless flavour 
tensors. The flavour tensors relevant to Ca — > e^7, are Sl,r and Tl^r, which respectively enccde flavour change among 
leptons of the same SU(2) representation, or chirality-changing flavour change. Formulae for these flavour tensors are 
in eqns (7) to (13). They involve dimensionless k and y functions of the mass ratio in the loop, which are given in 
the Appendix. Assuming that the new particle masses arc at the scale Anrfc > v„ these k and y functions can be 
expanded in v"^ /Aj^ppQ, because one of the masses is a Standard Modele mass. The coefiicients in this expansion will 
be combinations of coupling constants. The terms in this expansion can be identified to an Effective Field Theory 
[42, 43] expansion in I/Atvpfc- of operators which can contribute to the dipole interaction. These include the usual 
dipole operator ( eqn (14) and first term below) at 0(1/A^pp,^), plus gauge invariant operators with an arbitrary 
number of extra Higgs fields: 

(-.(6) _ (~,{8) 

.^^{£0H)a^''PRe^F^, + -^f^H^^ H{IpH)a^'' Prb^F^, + ... (44) 

^^NPFC ^^NPFC 

In this way of thinking, C'^-* is our invariant, and it will be a good approximation to the dipole vertex functions if 
higher order terms in eqn (44) can be neglected, as would be natural for Ajvp v. However, the New Physics scale 
should be ^ 10 TeV to see r e^7 at a Super-B Factory. 

Normally, for such "TeV-scale" New Physics, the coefficient C(®) is matched to the dipole vertex function at the scale where the new 
particles are removed from the theory. 
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Two models were considered, to explore the possibility that the coefficient matrix C*^^) is an invariant, and an 
adequate approximation to the dipole vertex function. In the case of a non-SUSY, TeV-scale type I seesaw, for which 
the Lagrangian is given in eqn (16), the new physics interacts with the doublet leptons, so contributes to Sl. Two 
types of singlet scalar leptoquark were also considered, respectively coupled to down or up-type quarks. If a generation 
number is attributed to the leptoquarks, then the coupling constants at the vertices on either side of the loop are 
three index tensors, and the invariant is obtained by summing generation on both internal lines. At 0{l/Aj^ppQ), 
the flavour tensors Sl.r due to leptoquarks or the seesaw are proportional to the invariant 

(see eqns (21), (30), and (31)), where is the mass-squared matrix of the new particles, and Ai and A2 are the new 
couplings appearing in the left and right side of the blobs on the top row of figure 1 (In the models considered in this 
paper, Ai = A2). 

The leptoquark could also interact with singlet and doublet fermions via the coupling matrices and A^. The 
chirality flip neccessaxy for the dipole can arise on the internal quark line, allowing Sq to generate the flavour tensor 
T represented in the second row of flgure 1. This gives invariants of the form 

XlYuM-^xI, , XryIm-^xI 

(see eqn (32)). Prom these cases, we can extrapolate that the flavour structure of C^^^ is obtained by multiplying 

one SM Yukawa matrix with two new flavoured couplings matrices, and possibly an inverse mass-squared matrix for 
flavoured new particles (a dimensionless invariant can still be obtained by normalising by the smallest flavoured new 
particle mass^ -^1 )• The form of the invariant can be obtained by inspecting the relevant Feynman diagrams: vertices 
give a coupling constant, and heavy propagators give M~^. 

It remains the question of whether the vertex functions ^l,r can reasonably be replaced by the coefficient C^^-* 
of dimension six operators, which is proportional to an invariant. As discussed after eqn (44), there are corrections 
proportional to w^/A^p^^^, which is < 1 by assumption. The perturbative expansion is therefore well defined. 
However, in flavour physics, there can be small mixing angles and hierarchical couplings, so the relevant question 
is whether the 0{v'^ / Aj^ppf-.) terms we wish to drop are small compared to the invariant that was retained. The 
answer is model-dependent. In the case of the leptoquark So, the C(v^/A^pp^) terms are multiplied by a down-type 
Yukawa coupling squared, and it is reasonable to neglect corrections ^ rai/M— . However, in the seesaw model, the 

So 

0(v'^ / A?^ppq) corrections are ~ m'^/M^ (for M a TeV-scale singlet neutrino mass), and for the leptoquark 5*0, they 
are ~ m^/Mg^. Such corrections could contribute an observable r — )• epj rate for new particle masses ^ 2 TeV (see 
eqns (4) and (5)). 

There can also arise logarithms of flavoured mass matrices, as for example in eqn (32). These could induce errors 
in the reconstruction of a flavoured matrix Y, because the eigenvalues of Y and Y{1 + logF) are different, even 
though they are diagonalised by the same matrix. This is unfortunate. However, maybe logarithmic cirrections can 
be considered "small" in a flrst approximation. 

Another stumbling block, in approximating the dipole vertex function as an invariant multiplied by an unknown 
constant, is that several distinct invariants could contribute to C*^^^. This was the case for the So leptoquark, which 
could induce Sl.r, and (see eqns (26) to (29)). If new spurious are few, or if new particles are rare, this 
problem may not arise. However, as studied in section 4.3, the leptoquark So can induce all four flavour tensors 
with similar contributions to t — > e/37, while remaining consistent with phenomenological constraints. It could be 
interesting to study whether there is a unique invariant in models, such as the supersymmetric seesaw, with a single 
new spurion for the charged leptons (the neutrino Yukawa matrix), but several new particles. 

Finally, the invariant could be interpreted as a term from an MFV expansion, as these terms are also constructed 
by multiplying flavoured matrices from the Lagrangian ("spurious"). However, there are many terms in an MFV 
expansion, multiplied by arbitrary ^ 1 constants. The aim here is to obtain only one invariant, in the hope of 
reconstructing the spurious, so some justification for rejecting more complicated invariants is required. We suppose 
that sums over flavour indices arise from loops, so are accompagnied by a l/(167r^), which could suppress invariants 
containing more "spurious" . In addition, Yukawa matrices have a Higgs leg, which may be contracted in a loop, which 
naively gives an additional l/(167r^) suppression. Alternatively, if the Higgs leg remains as a vev, this brings a factor 
of v/Anp. In summary, considering only the simplest invariants may be justified, provided that its reasonable to 
neglect higher order terms in both the loop expansion, and the Effective Field Theory expansion in v/A^p. 
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5.2 Testing a hiercirchy? 

The current bound from — >■ 67 of 

ifci^ < 2.4 X 10-1^ 

is almost four orders of magnitude more restrictive than the planned sensitivity of SuperB Factories to r — >■ 67 
and T ^ fij (see table 1). Pessimists could interpret that the coefficients of Icpton flavour violating operators are 
suppressed, so r ^ 67 and t ^ fij will be beyond the reach of Super B Factories. This paper takes the opposite view 
that SuperB Factories will see r ^ 67 and/or r i^j. From a model-building perspective, this could arise because 
flavour change is larger in the; third generation. From a phenomenological perspective, one can only learn about New 
Physics by seeing it, so to address the question "what can the — ^ e/37 decays tell us about flavour structure in the 
lepton sector?" , we assume that r — >■ 67 and/or r — > /X7 are observed. Furthermore, the smallest angle in the CKM 
matrix is Vub ~ 4 x 10~^, so a 10~^ suppression due to a small mixing angle squared is "standard". 

In practise, we study the less ambitious question of whether a hierarchical flavour tensor makes falsifiable pre- 
dictions, where the definition of hierarchical is given in section 4.1. It is model- independent, and unrelated to the 
invariants discussed in the first part of the paper. We focus on a hierarchical structure, because it has sufficiently few 
parameters to be predictive, and because the Yukawa matrices are hierachical. A hierarchical flavour tensor which 
induces r ^ 67 or r — >^ /X7 at a Super-B Factory, must contain a small parameter to suppress fi ^ ej below the 
current bound. Section 4 presents the argument that this small parameter also suppresses a r — > 6^7 decay, implying 
that a hierarchical flavour tensor forbids the observation of all the modes r — > 61,7, r — >■ eRj, t — >■ /xz,7 and r — >■ ^irJ. 
This argument rests upon five assumptions, given in section 4.2. The third and fifth appear problematic. The third 
assumption is that new particle mass scale is ^ 600 GcV, so that the couplings which control the observed t e/37 
rate are "large" . This enforces a hierarchy with respect to the couplings which control fj, — >■ 67, and the other r — >■ e'^j 
decay. The LHC may verify this assumption in several models. The fifth assumption is that the dipole vertex function 
should not have its general form given in eqn (3), but rather should be dominated by either S^ ji or T. This means 
that New Physics should be in one of the doublet or singlet lepton sectors (generating Sl.r)j or, if it interacts with 
both, then the amplitudes with chirality fiip inside the loop (T) should dominate those with chirality flip on the 
decaying fcrmion line {S^ ji). This condition arises because, for a given observed t — > e^7 process, S^.^f and T predict 
different rates to be suppressed, as is summarised in table 3. An example where this condition is not satisfled is given 
in section 4.3.2. However, if new particles below 10 TeV are rare, this flfth assumption is not unreasonable. 

Finally, if a hierarchy were measured, then it could be interesting to suppose that the flavour tensor is an invariant, 
and explore what this implies about the flavoured matrices of the Lagrangian. Section 4.3.1 briefly discusses the simpler 
question of when a hierarchy in the Lagrangian flavour structures is transmitted to the flavour tensors contributing to 
the dipole vertex function. 

5.3 Summary 

Beyond the Standard Model physics is required in the lepton sector for neutrino masses. If present below 10 TeV, it 
could induce observable e„ e/37 decays in upcoming experiments. Observing such decays could give information on 
the flavour structure of this New Physics. 

A hierarchical (in flavour space) dipole vertex function (where hierarchical means dominated by its largest eigen- 
value), could be confirmed if some, but not all of the r e^7 decays are observed. Were it possible to distinguish 
internal from external line chirality flip, such a hierarchy could be tested. 

In simple models (with few loop diagrams containing one SM particle), the coeSicient of the 0{1/M^gj^) term in 
the dipole vertex function is a "Jarlskog-like invariant" , meaning that it is obtained by multiplying flavoured matrices 
from the Lagrangian. It has the form MM^gj^^X^, where Ai and A2 are the flavoured couplings on the two sides of 
the loop. The invariant is a better approximation to the vertex function when the particles in the loop are (very) 
weakly coupled to the Higgs. Invariants elegantly give a linear relation to flavoured matrices in the Lagrangian, avoid 
confusion about basis, and may simplify the combination of bounds from various processes on New Physics. 
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Appendix: Summary of results from Lavoura [8] 

Lavoura defines coefficients aL,<JR, related as 
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to tlie Xx used here. The ax are written in terms of functions ki, kf, etc, which are m^x the definitions given here 
( nia is already scaled out of the definition of Sl.r in eqn (3)): 
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where * = rnj/m\, and niB is the boson mass. 
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